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Abstract 



The importance of the spectral density of the Dirac operator in study- 
ing spontaneous chiral symmetry breaking and anomalous U(1)a symmetry 
breaking are reviewed. It is shown that both types of symmetry breaking 
can be traced to effects of modes near zero virtuality. Above T c , where chiral 
symmetry is restored, it is shown on general grounds that (in the m q — > 
limit), the density of states vanishes at zero virtuality faster than A, where A 
is the virtuality — i.e. p(X) ~ |A| a is not possible for a < 1. Isospin invari- 
ance is used to show that p(X) ~ m*~ a |A| a is also not possible for a < 1. 
State-of-the-art lattice calculations are reviewed in light of these constraints. 
In particular, it is argued that violations of these constraints by lattice cal- 
culations indicate possible large systematic errors; this raises questions about 
U(1)a violating effects seen on the lattice. It is also shown that above T c , the 
Dirac spectrum has a gap near zero (in the m q — > limit) unless contributions 
from quark-line-connected and disconnected contributions conspire to cancel. 



I. INTRODUCTION 



It is a great pleasure to give this talk at a workshop in honor of Mannque Rho's 60 
birthday. Mannque has been a central figure in our field for quite a long time. When I 
began working in nuclear physics about 15 years ago, it became quite clear to me that 
Mannque was quite special. I soon realized that our field could be roughly divided into 
two classes — theoretical physicists with powerful and elegant mathematical tools at their 
disposal, and those with a phenomenological bent who had the intuition to make the kinds of 
simplifications needed to deal with the complications of realistic strongly interacting systems 
and make experimentally relevent predictions. Mannque is almost unique in belonging to 
both groups. 

Having divided the field this way, let me begin this talk with an apology. Virtually noth- 
ing which I will deal with here has immediate phenomenological consequences. The problem 
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with which I am dealing concerns the QCD vacuum in thermal equilibrium above the chiral 
phase transition. This raises the obvious question of just how does one heat up vacuum? 
An idealized experimental setup — as designed by a typical American student — consists of 
an empty test tube held above a large candle. I have been told that this method will not 
work in practice, and that the only practical way to explore this phase in the laboratory 
is through ultrarelativistic heavy ion collisions. Unfortunately, the way the physics of the 
phase of interest gets translated into experimental signatures which can be observed in heavy 
ion reactions is by no means clear. I will barely touch on this question in this talk, and what 
I do suggest will turn out to be highly sensitive to assumptions about the dynamics of the 
phase transition and will be highly speculative. 

A theorist might wish to attack a problem so far removed from experiment for reasons 
other than the sheer delight in being unable to be shown wrong by our grungy experimental 
friends. The fact is that we believe we know the underlying theory (QCD) and we can 
cleanly formulate well-posed questions about QCD in thermal equilibrium in terms of the 
degrees of freedom of QCD. Moreover, at least in principle, these questions can be answered 
to arbitrary precision via sufficiently good lattice simulations. Of course in practice, present 
day lattice simulations have some serious shortcomings for technical reasons associated with 
a lack of sufficient computer power, but there is every reason to believe that as computers 
get better, the lattice community will slowly converge on the right answers. 

Since we are ultimately concerned with issues of nonperturbative QCD, but we will be 
working with the QCD degrees of freedom — quarks and gluon-it is probably useful to intro- 
duce a few tools. One essential tool is the composite operator or interpolating field. These 
beasts are constructed from the quark and gluon operators. Typically one is interested in 
local gauge-invariant operators with fixed transformation properties under Lorentz trans- 
formations, parity, and isospin rotations. I will generically note such operators as J and 
examples include: 



which have the quantum numbers of a pseudoscalar isovector and scalar isoscalar, respec- 
tively. One can construct an infinite number of such operators. Essentially all of the infor- 
mation one can extract from QCD is encoded in the various n-point correlation functions 
of these currents. A typical example is the two-point function (J(x,t)J(0,0)) where one 
can choose the appropriate boundary conditions of interest (e.g., time-ordered, advanced, 
retarded, etc.) These correlators tell us about the response of the QCD vacuum given a 
disturbance associated with J. 

How can one calculate these things? In principle, there is a straightforward way to 
calculate the expectation value of any observable given in terms of quark and gluon field 
operators. The method is functional integration. Here we will discuss the problem for 
imaginary time (Euclidean space) and "all" one needs to do to get real time results is 
analytically continue. The Euclidean functional integral for the expectation value of any 
operator is given by 



J ~ qi^r q 




(1) 
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Z = J [DA] [Dq] [Dq] e i5 « CD (2) 

where Sqcd is the action. 

Here we are interested in QCD at finite temperatures. One of the virtues of the Euclidean 
space functional integral is that one can calculate finite thermal expectation values exactly 
as one calculates zero temperature ones. Indeed the functional integral expression in eqs. (0) 
holds for thermal expectation values. The only difference from the free-space case concerns 
the class of configurations over which one integrates. To use functional integrals at finite 
temperatures one simply imposes the boundary condition |J that all boson fields (gluons 
in this case) are periodic in time with a period given by /?, the inverse temperature, while 
fermion field configurations are anti-periodic: 

Apfat) = A„(x,t + P) q(x,t) = -q(x,t + /3) 

The above formulation holds for any operator. The operator which I will stress in this 
talk is the density of states of the Dirac operator; the question I will address is how this 
relates to various correlation functions. The operator can be defined in the following way: 

p(X) A =l/Vol £5(A-A,) (3) 

j 

where i\j is an eigenvector of the (Euclidean) covariant Dirac operator acting on spinors, 

D^Qj = i^jClj , 

and Vol is the volume of space-time with a finite temporal extent given by (3 and, ultimately, 
infinite spatial extent. Note that the Dirac operator D^j 11 is anti-Hermitian in Euclidean 
space and hence has purely imaginary eigenvalues. Note also that pa is implicitly a functional 
of the gauge configuration A through the covariant derivative. 

A useful property of the Dirac operator is that it anti-commutes with 75 which, in turn, 
implies that for every nonzero eigenvector qj, is an eigenvector with opposite eigenvalue. 
This, in turn, implies that Pa(A) is an even function of A. The parameter A will be referred 
to as the virtuality. There is an extremely simple way to represent pa as a trace which we 
will find useful later: 



Pa (A) = lim7r/Vol Im ( Tr 



A + iD^j^ — ie 



(4) 



where Tr indicates a trace over Dirac indices and also a functional trace over configurations. 

The thermal expectation of pa will be of prime interest. I will define p(A) to be expec- 
tation value: 



p(\) = (p(X)a) (5) 

It is not easy to gain an immediate intuitive grasp of the meaning of p(A). However, as will 
be stressed in this talk, p(A) in the vicinity of A = is intimately related to the spontaneous 
breaking of chiral symmetry and anomalous breaking of U(l) axial symmetry. 
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Before proceeding, let me remind you of some well-known facts about chiral symmetry 
in QCD. The up and down quark masses in QCD (which I will generically call m q ), are 
much less than the typical hadronic scales. Neglecting the quark masses is a reasonable 
first approximation which we will take here and the effects of nonzero quark masses can 
be computed via chiral perturbation theory. Now if we neglect quark masses the QCD 
Lagrangian is invariant under 

q -> e *»^/ a 

where f are the Pauli matrices and 9 are arbitrary parameters. While the underlying theory 
has this symmetry, the ground state (i.e. the vacuum) does not, at least for T < T c : the 
symmetry is spontaneously broken. One easy way to see this is to note that the chiral 
condensate (qq) is nonzero. This indicates spontaneously chiral symmetry breaking since 
under a chiral rotation through tc, in the a direction, the operator qq transforms into qi^T a q 
which has a zero expectation value. The consequence of this spontaneous breaking are 
profound: through Goldstone's theorem we see that the symmetry breaking implies the 
existence of massless isovector pseudoscalars which we interpret as pions. (In nature, they 
are merely very light and not massless since the quark masses are not exactly zero). Similarly, 
the spontaneous symmetry breaking implies that the pions don't interact at threshold (in 
the zero quark mass limit) which allows one to formulate a chiral perturbation theory for 
low energy interactions. J|] 

It is generally believed that in the high temperature phase of QCD that chiral symmetry 
is restored. That is to say that the symmetry of the thermal state is the same as the 
symmetry of the Lagrangian. One consequence of this is above T c the chiral condensate 
vanishes. 

Let me now turn to the connection between the chiral condensate and the spectrum of 
the Dirac operator. Banks and Casher proved a remarkable result in 1980: || 

lim (qq) = tt p(0) (6) 

m q — >0 

The chiral condensate in the zero quark mass limit is directly proportional to the density of 
states at zero virtuality. The proof of this is quite straightforward: 

(qq(x = 0)) = (tr[SU(0, 0)]} (7) 

where Sa is the quark propagator in the presence of a background gluon field, 

(D^ + m q )S A (x,y)=5 4 (x-y), (8) 

and tr indicates a trace over the Dirac indices only. Since Sa(x, y) is simply the inverse of 
the operator + m q , it is clear that spatially averaging over all x will simply yield the 

functional trace of (D^v- + m, 



-i. 



qj 



+ m 



d xtr[SA(x, x)] = Tr 
One can exploit translational invariance to note the 

r i 

) 



(<?<?) = / d 4 xqq(Q )) = l/Yol(7V 
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Taking the limit asm->0 and comparing with eq. (|4]) immediately yields the Banks- Casher 
relation of eq. @. 

One obvious and beautiful thing about the derivation is that it only depends on the 
Euclidean functional integral formulation of the problem and hence is automatically valid at 
finite temperature. Above the phase transition where chiral symmetry is restored and (qq) 
goes to zero (or more strictly 0{m q ) ) we see that p(0) — > (or more precisely 0{m q )). So 
we see that p(0) serves as an order parameter for chiral symmetry and its restoration. 

The spectral density of the Dirac operator also gives us useful insights into the U(1)a 
problem (i.e. the axial U(l)) problem). It is well known that the QCD Lagrangian with 
massless quarks has an additional axial symmetry: 

q _> e *»»/2 g (9) 

This symmetry is anomalously broken by quantum effects — the current associated with it 
is not conserved. The reason for this is quite simple: divergences in the theory require 
regularization; there is no way to regularize which simultaneously respects U(1)a symmetry 
and gauge symmetry. Since the renormalization of the theory requires unbroken gauge 
symmetry, U(1)a must break. The consequences of this are far reaching. For example, this 
is an essential ingredient in the understanding of why there is no U(1)a Goldstone boson. Q 
Now, just as there is a connection between the spectral density at zero virtuality and the 
spontaneous breaking of ordinary chiral symmetry, there is an important connection between 
the states at zero virtuality and the spontaneous/anomalous breaking of the U(1)a- The 
point is very simple — all effects associated with U(1)a breaking in massless QCD can be 
traced to modes of the Dirac operator in the neighborhood of zero virtuality. || To make 
this concrete, consider what I will call as the U{1)a susceptibility, Xu(i) A > which is simply the 
spatial integral of the correlator in the isovector pseudoscalar (jr) channel minus correlator 
in the isoscalar pseudoscalar ( " 77') " channel: 

Xu(i) A = J d 4 x (j n (x)j n (0) - jr,'{x)j v >{0)) (10) 

with j n = Tpi^rq and = q^q. An analysis very similar to that of the Banks-Casher 
relation immediately yields; 

2 

/771 

^From the form of eq. ( O ) it is obvious that in the limit m q —>■ the susceptibility must 
vanish unless there is strength in the immediate neighborhood of zero. The important thing 
to realize is that as was pointed out a couple of years back || this is generic: all U(1)a 
violating amplitudes get their strength entirely from the A ~ region. 

A few years ago, motivated by an insight of the previous sort (in the context of the 
instanton liquid model), Shuryak M asked the following provocative question: since we know 
above T c , p(0) = and we know that U{1)a violating amplitudes come from modes at A = 0, 
is it possible that all effects of the anomalous U(1)a breaking vanish in the chirally restored 
phase? At first sight this seems completely nuts — the anomaly is an operator equation and 
the axial current is not conserved regardless of the state of the system. On the other hand, 
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it is well known that it is only due to subtle interplay between topology and the anomaly 
that any effects of U(1)a breaking are seen. In fact, a couple of years ago I showed on very 
general grounds based on the positivity of the measure in the QCD along the lines of the 
QCD inequalities of Weingarten, Witten and Vafa @ that unless there are contributions to 
the functional integral which form a set of measure zero in the massless limit of the theory, 
U(1)a violating amplitudes must vanish above T c . This raises the obvious question as to 
whether such a set of measure zero can contribute. As noted by Lee and Hatsuda || and 
by Evans, Hsu and Schwetz 0, such contributions are possible in a dilute instanton gas 
where one might expect p(X) ~ m u m^5(A). Such a form clearly is a set of measure zero in 
the massless quark limit but would contribute to U(1)a violating amplitude. This leaves us 
with the question of whether QCD actually does yield U(1)a violating effects above T c or 
not. 

Now the preceding analysis is highly formal and mathematical and it is useful to remind 
ourselves of what Bertrand Russell once said: "Mathematics may be defined as the subject 
in which we never know what we are talking about." 

Given the bleak prospects of using purely formal reasoning to answer this, it seems there 
are two possible avenues of attack. The first is through phenomenology and the second is 
via lattice studies. 

As I noted at the outset it is very difficult to use phenomenology to pin things down 
since so much of the dynamics is unknown. There is at least one scenario where we can get 
a dramatic effect but it depends on many unproven assumptions. The first is that U(1)a 
violating effects do vanish above T c ; this is, of course, what we are trying to check. The 
second is that the transition is second order and happens sufficiently rapidly to get out 
of thermal equilibrium over large domains and develops an instability for pionic growth. 
This is the so-called disoriented chiral condensate scheme of Wilczek and Rajagopal [JTO 



If this happens, then along with collective enhancements of low Pt pions we would get 
enhancements of low Pt rj's. Clearly this scenario is highly speculative and we turn to the 
question of what the lattice can tell us. 



II. LATTICE RESULTS 

In principle, the lattice can answer our question about U(1)a violations above T c . Of 
course, there may be the usual problems of interpretations given effects of finite masses, 
lattice spacings and lattice sizes and the like. While there have been a number of lattice 
calculations of quantities sensitive to U(1)a violations many of them depend on very nu- 
merically unstable quantities such as screening masses in various channels which we know 
are highly sensitive to quark mass effects and threshold effects. A more natural quantity to 
study is the U(1)a violating susceptibility defined in eq. ([0]). To my knowledge there are 
only two calculations in the literature of this quantity. The first by Chandrasekharan and 



Christ [|TTJ and the second by Bernard et al. [JL2] . Both of these papers claim that that they 
see evidence that Xu(i) A does not go to zero above T c indicating that U{1)a violating effects 
do not vanish above T c . In particular, Chandrasekharan and Christ claim that they see a 
small but nonvanishing Xu(i) A f° r modestly small quark masses. Bernard et al. did a more 
systematic study; they calculate for a number of modestly small quark masses and extrap- 
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olate back to zero and obtain a nonzero value. In contrast, the ordinary SU(2) x SU(2) 
chiral symmetry breaking susceptibility (the pion-sigma channels) goes nicely to zero. 

To some this may seem to settle the question; but to quote Mark Twain, 11 It's differences 
of opinion which make horse races." There are several reasons to doubt the lattice studies. 
The first is the obvious one; the calculation involves a small number from an extrapolation of 
quantities with unknown and potentially large systematic errors due to finite masses, lattice 
spacings sizes etc. Moreover at least one class of these systematic errors (finite quark masses 
compounded by finite lattice spacing) can easily account for a spurious U(1)a violating 
amplitude without affecting ordinary chiral symmetry. The reason is quite simple: due to 
the fermion doubling problem, viable lattice schemes in general violate axial symmetries in 
a spurious way. There are different schemes for dealing with fermion doubling which all 
have the general feature that the spurious axial symmetry violating effects vanish in the 
continuum limit of a — > 0. The calculations by Chandrasekharan and Christ and Bernard et 
al. used staggered fermions. These have the property that even at finite lattice spacing one 
axial symmetry is conserved (for m q = 0). This conserved axial current is the one which is 
conventionally associated with the pion channel. On the other hand, the U(1)a symmetry is 
explicitly broken by the lattice formulation and is only restored in the continuum limit. Thus, 
for example, even in the case of noninteracting massless fermions in this lattice formulation 
one would see U{1)a violations until a is sent to zero. 

Accordingly it would be very useful to study the anatomy of these calculations to see 
whether the results are coming from lattice artifacts. As was discussed in the introduction, 
Xu{i) A i s completely determined by p(A) near zero. Thus if one could measure p(A) for the 
lattice configurations contributing, one could see from where the U(1)a violating amplitudes 
get their strength. In fact, as we will see, the U(1)a violations above T c in ref. 11 get their 
strength in a manner which appears to be inconsistent with chiral restoration suggesting the 
result may well be due to lattice artifacts. 

While a direct lattice calculation of p(X) is what we want, there are presently no available 
such calculations. This is not surprising as it would be a numerical nightmare to calculate. 
Fortunately an integral transform of p(A) has been calculated by Chandrasekharan and 
Christ. (Tl| In particular they calculate 



^> = / dA '< A >«T^ (12) 

This is actually a relatively easy calculation to do on the lattice since it amounts to calcu- 
lating (qq) with a quark mass of in the propagator instead of m q . Hence standard codes 
can be used. Since the quark mass in the propagator can be made much less than that in the 
functional determinant (which is implicitly contained in p) one can call such a calculation 
partially quenched. Note that by construction, f{m q ) = (qq). 

Chandrasekharan and Christ calculated f{m^) for various choices of coupling constants 
which correspond to different values of the temperature. (Let me remind you that by 
changing the coupling constant with fixed physical observables one is effectively changing 
the lattice spacing. Since the number of steps in the temporal direction is fixed, this in turn 
changes the periodicity of the lattice in time and hence the temperature.) The important 
thing for our purposes is that for temperatures above T c and small m^, the function is 
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apparently linear on a log- log plot over several decades implying that f{m^) ~ m^ for 
sufficiently small A. Moreover it is found that in all cases a < 1. 
Now it is easy to see that this behavior implies 

p(\) = c\X\ a , 

where c is a constant, and this form applies at sufficient small A. The simplest way to see this 
is to put this form into eq. (0) and evaluate the integral. One finds that p ~ |A| a implying 
/ ~ m^. Indeed this is seen easily by dimensional analysis. In the following section, I 
will show that there is a constraint imposed by chiral restoration which shows that a form 
p ~ |A| Q with a < 1 cannot occur. This in turn raises questions as to whether the lattice 
results are dominated by artifacts. 



III. CONSTRAINTS 

Recall that we are studying the symmetry restored phase. In a chiral symmetric phase 
one knows that the log of the partition function is an even analytic function of the quark 
mass near zero. The reason for this is quite easy to understand. Look at the n th derivative 
of log(Z) with respect to the quark masses: 



d n log(Z) 



dm q 



=o V 



1 (l(d 4 xqqr). (13) 



In the chiral restored phase this quantity should not change under a chiral rotation. Under 
a chiral rotation through it, qq — > — qq, thus the preceding expression must vanish for all 
odd n; in general, for even n we expect nonzero finite values implying that log(Z) is an even 
analytic function of m q . 

We should recall that the quark condensate is given by 

i aiog(z) 

( ^ = voT^T' (14) 

and that from the Banks-Casher analysis, 
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Combining these we see that 



h'^M^/^JiL (16) 

m,-o Vol dm n q dm™- 1 J A 2 + m 2 



Recall that this must be zero for all odd n. 

Before proceeding with a careful analysis, let me begin with a small swindle. For your 
convenience, I will identify the swindle at the outset: assume that in eq. ([16]) the quark mass 
derivative acts only on m q /X 2 + m 2 , and not on p(A). (I should remind you that p depends 
implicitly on m q through the functional determinant.) Diagrammatically, this corresponds 
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to including only quark-line-connected diagrams. I will label any quantity calculated at this 
level with the superscript qlc. From eq. flT6| ) we see that 

d n log(Z) 

hm 



+ o dm™ 



lim / dX p(X) r - q = = - / dXp(X) — . (17) 

-»,->o J FV ' dmT 1 X 2 + ml J w A n y ' 



q q 



Now if p(A) ~ |A|| a for small A, then from eq. (|T7|) one sees that the integral diverges for 
any n > a. On the other hand if log(Z) is an even analytic function this is impossible. We 
conclude from this analysis the following: if only quark-line-connected diagrams contribute, 
then p (for m q = 0) is infinitely flat in A — it goes to zero faster than any power law. This is 
what you would expect if p had a gap at zero. By a gap at zero, I simply mean that p(A) = 
unless A is greater than some minimum. Alternatively, one could imagine a situation in which 
p had an essential singularity at A = 0, such as p ~ e _1 /\ In any case, p is infinitely flat at 
the origin. 

My guess is that the preceding description is correct — that above T c , p does, in fact, 
have a gap. I cannot prove it however. The problem is that the preceding analysis only 
included the quark-line-connected diagrams. It is possible in principle that the quark-line- 
disconnected parts, coming from derivatives of m q acting on p(A), conspire to cancel the 

divergent quark-line-connected part. One can prove, however, in the chiral restored phase, 

dp 



that for m q = 0, p must go to zero faster than linearly: . A 



= 0. 

A=0 

The proof of this constraint is actually quite simple. Consider the susceptibility associ- 
ated with the pion channel, 

X*= I d A x (j n (x)j n (0)) , (18) 



with j w = q^rq. By standard analysis of the propagator analogous to the derivation of the 
Banks-Casher result discussed in the introduction, one finds: 

x. = / = M < 19 > 

J X A + m l q m q 

This last form is just a functional integral derivation of the Ward identity implicit in the 
Gell-Mann-Oakes-Renner relation. Now suppose that for small A, p = c\X\ a with a < 1, and 
c is a constant of proportionality. It is a trivial exercise to evaluate the integral in definition 
of Xn using this form, and one finds that 

X, ~ cm a -\ (20) 

Let us suppose for the moment that the form p = c\X\ a survives in the chiral limit of m q — > 
so that c is finite in this limit. Evaluating the integral in eq. ( p0|) for this case gives for 
small m q 

7T 

Xn = J J— (21) 

m v q a cos(a7r/2) 

which diverges in the chiral limit m q — > 0, provided a < 1. (For a = 1 it diverges as 
log(m g ).) Moreover, if as hypothesized, one were studying the chiral restored phase then as 
m q — > the susceptibilities in the pion and a channels must be equal. On the other hand 
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Thus, if Xtt (and hence Xa) diverge in the chiral limit we see that log(Z) is not an analytic 
function of m q in contradiction to the system being in a chiral restored phase. 

^From the preceding analysis one can deduce that in the chiral limit of m q — > and the 
chirally restored phase, p(A) cannot go as c|A| a for a < 1. This appears to be inconsistent 
with the lattice results of Chandrasekharan and Christ. Of course, the preceding analysis 
was done for the chiral limit of m q — > 0, while the lattice studies were done with a fixed 
quark mass. This raises the following possibility. Suppose the constant c used above were 
not a constant but depended on the quark mass in such a way that it went to zero as 
m q — > 0. Would such a scenario allow a U(1)a violating susceptibility while still being 
constant with the observation that f{m^) ~ m^? At first sight it appears to be easy to do 
this: /(£) ~ m q ~ a m^ 1 i.e. the constant c could be proportional to m q . This possibility is 
consistent with the lattice data in Chandrasekharan and Christ. Since they do not do the 
calculation for multiple values of m q at fixed it is impossible to tell if this behavior is 
present. If f{m^) behaves this way it is clear that for small A and m q , 

p(A) = bm\~ a \\\ a (23) 

with b a constant. If p is given by eq. ([23|) it is a simple matter to verify that Xtt an d Xa are 
finite and equal and that Xu(i) A — (1 — a )Xw 

The preceding scenario appears to show a way to reconcile the lattice data with chiral 
restoration and U(1)a violation. In fact, this scenario is inconsistent, unless there is spon- 
taneous symmetry breaking of isospin and is therefore not physically realized. This leaves 
an inconsistency which suggests that the lattice results may be dominated by unphysical 
lattice artifacts. The issue of isospin has up until now been ignored. Implicitly I have taken 
the mass of the two light quarks to be degenerate. Assuming isospin is not spontaneously 
broken, this is not an issue since at the end of the problem both the up and down quark 
masses are taken to zero. If isospin is spontaneously broken, the order in which the limit is 
taken can matter, but if it is not spontaneously broken, then the ordering is irrelevant and 
we might just as well take the masses equal. The problem with the form in eq. ( p3"D is that 
it depends on the quantity m q ~ a . If one generalizes to m u ^ one must generalize m q ~ a 
into some function of m u and irtd which goes to m]^ a when m u = rrid = m q . Consider the 
following quantity: 



Xud = J d A x (uu(x) dd(0)) 
Vol dm u drrid 

Plugging in p from eq. fl23|) into a generalization of eq. fll6|), for the case where m u and 
are varied independently, and evaluating the integral yields: 

dm q (m u ,m d ) 1 - a 6vrm° 

Xud — o 7 77TT • l^J 

orrid cos(a7r/2) 
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If isospin and chiral symmetry are both unbroken then the result for \ ud will be finite as 
and is independent of how the quark masses go to zero. Let us assume this to 
be true and show that it leads to a contradiction. On dimensional grounds we can always 
write 

—1— = (1 - a)m~ alpha g(m u /m d ) (26) 
om d 

where g is some presently unknown function. Inserting this form into eq.(E5|) gives 



Xud = g(m u /m d ) — j— (27) 

cos(a7r/2) 

By hypothesis, isospin and chiral symmetry are not spontaneously broken so that the value 
of Xud must be the same independently of how the chiral limit is approached. This in 
turn means that the value of Xu,d has to be independent of the ratio m u jm d since we can 
approach the chiral limit with this ratio fixed to any value we wish. From this we deduce 
that consistency requires g = 1 and m}~ a = m^rrQa which in turn implies 

p(A) = bm d m-a \X\ a (28) 

At this point, however, we see a contradiction: the preceding expression is not invariant 
under m u <-> m d ; on the other hand since theory is isospin symmetric up to the values of 
the quark masses the spectral density must be invariant under the switching of the up and 
down quark masses. Thus we conclude that the scenario in which p(X) ~ m q ~ a \X\ a is not 
consistent. 



IV. SUMMARY 

Let me summarize the situation: I have discussed the concept of the spectral density near 
A = 0, and argued that it was interesting and an important theoretical tool for understanding 
both spontaneous chiral symmetry breaking and the spontaneous/anomolous breaking of 
U(1)a- One problem of physical interest that this tool may prove useful for is the question 
of whether the effects of U(1)a- axial symmetry breaking are manifest above the chiral 
phase transition. Although present lattice simulations seem to indicate that U(1)a violating 
effects survive above the transition, it is possible that the calculations are dominated by 
lattice artifacts. One way to see if this is so, is to test whether the spectral densities implied 
by the lattice calculation are consistent with constraints imposed on the spectral density by 
chiral symmetry and isospin in the unbroken phase. 

There are a number of such constraints: first, one sees that unless quark-line- 
disconnected graphs conspire to cancel the quark-line-connected ones, in the chiral limit, 
p(X) is infinitely flat in the sense of having all derivatives vanish as one would expect from 
a spectrum with a gap. Second, in the chiral limit, || = 0, indicating a behavior such as 
p(X, m q = 0) ~ |A| a , is not possible for a < 1. Finally it is shown that for small A and m q a 
behavior of the form p(X,m q = 0) ~ m} q ~ a \X\ a is not possible for a < 1. 

Given these constraints it is reasonable to ask whether the lattice calculation of Chan- 
drasekharan and Christ, which have p behaving like a power law with a < 1 (and hence 
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in apparent violation of the constraint) are dominated by lattice artifacts. Since the lattice 
calculations of ref. 12 do not study the spectral density near zero it is impossible to know 
whether they are consistent with the constraints derived above. However, as these calcula- 
tions were qualitatively similar to those in ref. 11 it is reasonable to question whether the 
spectral functions are consistent with being in a chirally restored phase. I conclude with a 
word of warning to potential consumers of lattice calculations for these problems: Let the 
buyer beware. 
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